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ABSTRACT 

In the pre-reionization Universe, the regions of the inter-galactic medium (IGM) which 
are far from luminous sources are the last to undergo reionization. Until then, they 
should be scarcely affected by stellar radiation; instead, the X-ray emission from an 
early black hole (BH) population can have much larger influence. We investigate the 
effects of such emission, looking at a number of BH model populations (differing for 
the cosmological density evolution of BHs, the BH properties, and the spectral energy 
distribution of the BH emission) . We find that BH radiation can easily heat the IGM 
to 10'^ — 10* K, while achieving partial ionization. The most interesting consequence of 
this heating is that BHs are expected to induce a 21-cm signal (JTf, '^ 20 — 30 mK at 
z ^ 12) which should be observable with forthcoming experiments (e.g. LOFAR). We 
also find that at z ;S 10 BH emission strongly increases the critical mass separating 
star-forming and non-star-forming halos. 

Key words: diffuse radiation; cosmology: theory; black hole physics; galaxies: for- 
mation; intergalactic medium 



1 INTRODUCTION 

The Sloan Digital Sky Survey^ (SDSS) has unveiled the 
existence of quasars at redshift z ^ 6 (Becker et al. 2001; 
Djorgovski et al. 2001; Fan et al. 2001, 2002, 2003; White 
et al. 2003). This indicates that super-massive black holes 
(SMBHs) with a mass of 10®"^° M© had already formed 
when the Universe was less than 1 Gyr old (Fan et al. 2001, 
2003). 

The processes which lead to the formation of such huge 
black holes (BHs) already in the early stages of the life of 
the Universe are very uncertain. A possible scenario is that 
SMBHs were built up starting from a seed intermediate- 
mass BH (IMBH, i.e. a BH with mass of 20-10^ Mq), which 
increased its mass by accreting gas and/or by merging with 
other IMBHs. 

In particular, if first stars are very massive (> 260 M©) 
their fate is to directly collapse into BHs, nearly without los- 
ing mass (Heger & Woosley 2002). This can produce a popu- 
lation of IMBHs, which are expected to efficiently accrete gas 
in the high-density primordial Universe and eventually to 
coalesce with other BHs (Volonteri, Haardt & Madau 2002, 
2003; Islam, Taylor & Silk 2003, 2004; Volonteri & Perna 
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2005; Volonteri & Rees 2005, 2006). Furthermore, the accre- 
tion of these IMBHs might be enhanced also during galaxy 
mergers, which tend to drive gas into the inner regions of 
the host galaxy (Madau et al. 2004). However, recent sim- 
ulations by Pelupessy, Di Matteo & Ciardi (2007) suggest 
that the accretion history of such seed IMBHs can hardly 
account for the SMBHs of the SDSS. 

On the other hand, seed BHs can be produced also 
by the direct collapse of dense, low angular momentum 
gas (Haehnelt & Rees 1993; Umemura, Loeb & Turner 
1993; Loeb & Rasio 1994; Eisenstein & Loeb 1995; Bromm 
& Loeb 2003), driven by turbulence (Eisenstein & Loeb 
1995) or gravitational instabilities (Koushiappas, Bullock 
& Dekel 2004; Begelman, Volonteri & Rees 2006, hereafter 
BVR06; Lodato & Natarajan 2006). In particular, the so- 
called 'bars within bars' mechanism (Shlosman, Frank & 
Begelman 1989; Shlosman, Frank & Begelman 1990) im- 
plies that bars, which form in self-gravitating clouds under 
some assumptions, can transport angular momentum out- 
wards on a dynamical time-scale via gravitational and hy- 
drodynamical torques, allowing the radius to shrink. This 
shrinking produces greater instability and the process cas- 
cades. BVR06 show that this process leads to the formation 
of a 'quasi-star', which rapidly collapses into a ~ 20 Mq BH 
at the center of the halo. The BH should encounter very 
rapid growth due to efficient gas accretion. 
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This allows the formation of seed BHs with mass 
< 10*^ Mq (BVR06; Lodato & Natarajan 2006, 2007), de- 
pending on the initial parameters (e.g. the temperature of 
the gas and the spin parameter of the parent halo). 

If such seed BHs formed at high rodsliift {z ^ 10 — 
30), they likely played a crucial role in the early Universe. 
Previous studies (Machachek, Bryan & Abel 2003; Madau et 
al. 2004; Ricotti & Ostriker 2004, hereafter RO04; Ricotti, 
Ostriker & Gnedin 2005, and references therein) showed that 
IMBHs accreting as miniquasars could be important sources 
of partial, early reionization. The efficiency of miniquasars 
in reionizing the high redshift Universe is especially due to 
the hardness of their spectra, which extend up to the X-ray 
band. 

For this reason, miniquasars are also indicated as 
sources of the X-ray background, and their density can be 
strongly constrained by the level of the unresolved fraction of 
this background (RO04; Dijkstra, Haiman & Locb 2004; Ri- 
cotti et al. 2005; Salvaterra, Haardt & Ferrara 2005; Volon- 
teri, Salvaterra & Haardt 2006; Salvaterra, Haardt & Volon- 
teri 2007). 

Finally, miniquasars can also heat the intergalactic 
medium (IGM; sec e.g. Nusscr 2005; Zaroubi et al. 2007, 
hereafter Z07+; Thomas & Zaroubi 2007), influencing a 
plethora of processes (the 21 cm line emission/absorption, 
the formation of the first structures, etc.). 

In this paper we analyze all the main effects that pri- 
mordial miniquasars can produce on the neutral IGM, i.e. 
on the regions of the Universe outside the ionized regions 
produced by the first BHs and stars. This is done by means 
of semi-analytical models combined with hydro-dynamical 
simulations. We consider all the most significant models 
for miniquasar formation, density evolution and spectra. 
WhereEis previous studies mostly focused on single aspects, 
our aim is to give a global description, as complete as pos- 
sible, of the role played by IMBHs in the early Universe. 

In particular, in Section 2 we will present an estimate 
of the radiation background produced by primordial BHs at 
high redshift. Then, we will discuss its effects on the IGM 
evolution, and in particular on the cosmic microwave back- 
ground (CMB) spectra and on the 21 cm emission (Section 
3). In Section 4 we will discuss whether the radiation back- 
ground can delay structure formation. In Section 5 we will 
discuss the relevance of our findings in the light of previous 
results. Finally, our results will be summarized in Section 6. 

We adopt the best-fitting cosmological parameters after 
the 3-yr WMAP results (Spergel et al. 2007), i.e. Ob = 0.042, 

= 0.24, Qdm = Om - Ob = 0.198, Oa = 0.76, h = 0.73, 
Ho = 100 /i km s"^ Mpc-^ 



ENERGY INJECTION INTO THE NEUTRAL 
IGM 



We start from considering that the mean free path of a 
photon of energy E emitted at redshift z is: 



X{E,z) = [nB{z)a{E)]- 



(1) 



where nsiz) = nB(0)(l -|- z)^ is the cosmological baryon 
number density at redshift z [fiB(O) ~ 2.5 x lO"'^ cm"''; 
Spergel et al. 2007], and a{E) is the photo- ionization cross 
section per baryon of the cosmological mixture of H and He, 
which is approximately 

( 0.75a h{E) 13.6eV < S < 25eV 

a(E) = l (T250 [-B/(250eV)]-2-®^ 25eV<S< 250eV (2) 
I ff25o[£^/(250eV)]-3-3 250eV<£;, 

where aniE) is the photo- ionization cross section of hy- 
drogen (see eq. 2.4 of Osterbrock 1989), and 0-250 — 3.2 x 
10~^^ cm^ is the cross section for 250 eV photons (see 
Zdziarski & Svensson 1989 for further details on the cross 
section at -E > 25 eV). In this paper we will neglect the 
absorption of photons with E < 13.6 eV. It is important 
to note that the above cross section is appropriate only for 
a neutral gas. Since the IGM close to luminous sources is 
mostly ionized (and the ionized fraction is not zero even in 
mostly neutral regions), eq. (1) might lead to an underes- 
timation of A, but this effect is important only in the last 
phases of reionization. 

On the other hand, the average distance between 'ac- 
tive' BHs is 



PBH(z)y 
C{Mbh){z) 



-1/3 



(3) 



where C accounts for the clustering of BHs^. 

The comparison of A and d shows that for photons of 
sufficiently high energy the mean free path can easily exceed 
d. For instance, the mean free path of a 500 eV photon 
emitted at a = 20 is about 9 comoving Mpc (~ d in typical 
models), but at the same redshift a 1 keV photon typically 
propagates for ~ 90 comoving Mpc. Since BHs are believed 
to emit a significant fraction of their luminosities at such 
energies, they will build up a roughly uniform background 
radiation field. 

The BH emissivity can be estimated by assuming that 
during active phases of accretion, each primordial BH pro- 
duces radiation at a fraction rj of the Eddington luminosity 
[LEdd ^ 1.3 X 10^^ ergs"^ (Mbh/Mq)], and that their av- 
erage spectrum at redshift z is described by some function 
F{E,z). Then, the proper emissivity is 



j{E,z)=C 



F{E,z) 



J F{E', z)dE 



-PBH,Q{z){jr^y{l + zf, (4) 



where C ~ 4.4 x 10 erg s ^ cm ^ Mg ^ is a normalization 



2.1 The total energy input 

First of all, we develop a formalism for estimating the total 
energy input of BHs into the neutral IGM at a given redshift, 
starting from basic properties of the BH population, taken 
from semi-analytic models (see next Section), such as the 
BH mass density psH at redshift z, the average BH mass 
(Mbh) at redshift z, and the duty cycle y of single BHs. 



^ If primordial BHs form in clusters (typically) of A^^; BHs, the 
probability that at least one of them is in an active state is 1 — 
(1 - (rather than y); so C ~ N^iy/ll - (1 - J/)^=']. In 

the following we will use C = 10 for models where BH clustering 
should be strong; this corresponds to N^i ~ 330(160) for y = 
0.03(0.06). 
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constant"^, and Pbh,q){z) is the BH density at redshift z, 
expressed in solar masses per cubic comoving Mpc. 

The mean specific intensity of the radiation background 
at the observed energy i?, as seen by an observer at redshift 
z, is then (cfr. eq. 2 of Haardt & Madau 1996) 

where the cosmological proper line element at redshift z' is 
dl c 1 1 



Table 1. Parameters of the BH growth model histories. 



dz' Ho 1 + z' [nuii + z'Y + nA]V2 ' 



(6) 



and r = t(E,z,z') is the optical depth effectively crossed 
by a photon emitted at redshift z' and reaching redshift z 
with an energy 



^^E,z,z') = jy~zf^.iE\±l 



)nB{z). 



(7) 



The definition of a cosmological background would re- 
quire that in eq. (5) Az — 0; but this is not appropriate 
for our purposes. In fact, as we already mentioned, we will 
be looking at regions outside the ionized 'bubbles' produced 
by the first luminous sources. So, we will examine the ef- 
fects of the radiation background on baryons quite removed 
from any particular source (we will refer to such baryons as 
the 'neutral-IGM' baryons), i.e. at a distance of the order 
of d/2. This is irrelevant for photons with long mean free 
paths {X{E, z) ^ d/2), but is of fundamental importance for 
photons with A <C d/2, which are absorbed in the vicinity of 
the BHs. In short, we start the redshift integration in eq. (5) 
from z+Az (where Az is chosen so as to skip the integration 
over distances < d/2), rather than from z. 

EYom the background spectrum 4>{E, z) we can easily 
obtain the energy input per baryon due to the absorption of 
background photons at redshift z, 



I 



e{z)=A-K I dEJ{E,z)a{E). 



(8) 



It must be noted that our use of the cross section (2) in 
the estimates of r (eq. 7) and of e (eq. 8) might induce two 
opposite errors. First of all, t is overestimated (and J under- 
estimated) when a significant fraction of the cosmic volume 
is ionized. On the other hand, when the absorbing medium 
is not completely neutral, we overestimate the fraction of ra^- 
diative energy which is actually intercepted by the baryons. 
The former effect leads to a moderate underestimation of 
e, starting at relatively high redshifts; the latter might lead 
to a large overestimation of e, but only for models where 
the IGM ionized fraction becomes quite large. We neglect 
both effects in our calculations: our results will generally be 
mild underestimates of the BH effects, except in the cases 
where the ionized fraction becomes large (a condition where 
we will significantly overestimate the BH effects). 

The energy input must be split into a fraction fio„ going 
into ionizations, a fraction fheat going into heating, and a 
fraction fexc going into excitations. These fractions actually 



^ C is chosen to be the emissivity (per cm^) when the BH den- 
sity is equal to IMq per comoving Mpc^ and BHs are assumed 
to accrete with efficiency 0.1. So, it is equal to 0.1 x 1.3 x 
10^8 (erg s-iM-i) x [3.086 x 1024(cm/Mpc)]-3. 



Model 


Pbh 


(Mbh)'' 


y" 


C 


Zs 


IMBH-3% 


106.75- 


-0.275Z 


103.5-0. 062 


0.03 


10 


30 


IMBH-6% 


107.6- 


-O.Sz 


106-0. 


0.06 


10 


30 


SMBH-3% 


107.26- 


-0.3752 


106. 5-0.062 


0.03 


1 


30 


BVR06'' 


106.675- 


-0.1876Z 


10« 


0.10 


1 


18 



° In solar masses per cubic comoving Mpc. 
' In solar masses. 

Consistent with the assumptions of the underlying models. 

Fig. 2 of BVR06 does not show the pbh evolution for 2 < 10; 
for this reason, at such redshift we will extrapolate the above 
formula. 



depend on the energy E of the absorbed photon; but ShuU & 
van Steenberg (1985) determined that, for all E ^ 100 eV, 
they are reasonably fitted by the expressions* 



fheat = 0.9971[1- 

fion = 0.3908(1 
f^xc = 0.4766(1 



0.2663n1.3163-] 

0.4092x1.7592 
Xh ) 
„0.2735\1.5221 



(9) 
(10) 

(11) 



where xh is the hydrogen ionization fraction {xh = 
n{H+)/[n{H°) + n{H+)]). As can be seen in Figs. 1 and 2, 
the contribution of photons with E < 100 eV to the back- 
ground is small or negligible (in the absence of reionization, 
the mean free path of 100 cV photons exceeds ~ 1 comoving 
Mpc and becomes comparable with d only at 2; ^ 5): so the 
use of these energy-independent functions is legitimate. 



2.2 Model parameters 

In the above section we have seen how we can obtain an 
estimate of the cosmological X-ray background produced by 
primordial BHs, and of the energy it can inject in the IGM. 
Such estimate mainly depends on three input quantities: the 
evolution of the cosmological density of BHs pbh{z), the 
duty cycle y, and the typical spectral shape of an active BH, 
F(E). The evolution of the BH average mass {Mbh){z), and 
the clustering factor C have much smaller effects. 



2.2.1 BH growth history 

There exist several models (e.g. RO04; BVR06; Z07-I-) pre- 
dicting the evolution of the BH mass density in the early 
Universe. Here we will discuss three different histories which 
are reasonable approximations of the models IMBH-3%, 
IMBH-6%, and SMBH-3% discussed in Z07-F (see their fig. 
8), and of one of the models in BVR06 (duty cycle 0.1, Mes- 
tel disk; from their fig. 2). The two IMBH models (IMBH- 
3% and IMBH-6%) assume that primordial BHs with mass 
~ IOOM0 form in small (10'^ - lO' Mq) and numerous halos, 
where H2 cooUng is efficient; the SMBH-3% and the BVR06 
models, instead, assume that primordial BHs of large mass 
( ^ IQ^Mq) form in larger (10*^ - lO^M©) and rarer halos 



* We report the expressions which are given for H, and neglect 
the small correction due to the presence of He. 
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cooled by atomic H. In all the four cases we will adopt the 
simple power-law approximations of the Z07+ and BVR06 
results which are given in Table 1. Such power-laws provide 
good fits to all the original models for z > 10, whereas at 
lower redshifts they are either a reasonable extrapolation 
(for BVR06), or give a slight-to- moderate underestimate of 
the predictions of the Z07+ models. Table 1 lists also the 
other parameters defining the BH growth histories: the duty 
cycle y is by far the most important, whereas our results are 
relatively insensitive to the assumptions on {Mbh){z) and 
C. This is quite fortunate as the value of y is intrinsic in our 
reference models, whereas none of them provides a simple es- 
timate of the other parameters. The values given in Table 1 
are guesses based on the general properties of the reference 
models, and on the notion that the typical BH mass should 
increase with time (especially when pbh grows fast). 

Furthermore, we assume that before a certain redshift 
Zs (zs — 30 for the evolutions taken from Z07+, Zs — 18 for 
the one from BVR06) the BH density (and emissivity) is 0. 



2.2.2 BH spectral energy distribution 

We experiment with three different types of BH spectral 
energy distributions (SEDs): simple power-laws F(E,z) = 
Fa{E), a template F{E,z) = Fsos,a{E) introduced by 
Sazonov, Ostriker & Sunyaev (2004), and a multi-component 
spectrum which is the sum of a multi-color black body and 
a power-law spectrum (see Shakura & Sunyaev 1973, and 
Salvaterra et al. 2005), F{E,z) = Fmc,-s-{E,z). 

Power-laws are characterized by their slope a, and are 
assumed to be 



<n 10 



F^{E) = 



0.01 eV <E < 10^ eV 
otherwise. 



(12) 



In the following, we will consider the power-law with a = 1 
(hereafter PLl SED) as our reference spectrum. 

The spectral template by Sazonov et al. (2004) is char- 
acterized by the slope in the 1 — 100 keV range, and its exact 
shape is 



FsosAE) 




0.1 eV < S < 10 eV 

10 eV < £ < 10^ eV 
10^ eV < _B < 10^ eV 
10^ eV < £; < lO'* eV 
otherwise, 



(13) 



where the constants Ci = iQ^i^-^-r) .^^^ ^ -^Q2.9-2a 
are chosen so as to ensure continuity, and the constant 
Co — 0.1607 ensures that the fraction of the BH luminosity 
which goes into photons with < 10 eV is the same as 
in the complete Sazonov et al. (2004) template (i.e. about 
0.85), even though we are not interested in the details of 
their model for i? < 10 eV. In this paper we considered the 
case a = 1 (SOSl SED): we chose such a relatively steep 
value (Sazonov et al. 2004 suggest values of about 0.7 — 0.8) 
because such SED is intended to show what happens with 
the steepest spectrum reasonably expected from BHs. How- 
ever, our results depend only weakly on this index. 

Finally, in the multi-component SED, the multi-color 
black body component is oc E^^^ and dominates up to a 
peak energy 




w 10 



10' 10' 10" 10= 10' 10' 10" 10= 
Photon energy [eV] Photon energy [eV] 



Figure 1. Spectrum of the background produced by primordial 
accreting BHs and seen by a neutral-IGM baryon in the four BH 
growth history scenarios (IMBH-3%: top left panel; IMBH-6%: 
top right panel; SMBH-3%: bottom left panel; BVR06: bottom 
right panel) at various redshifts (z = 8: thick dotted line; z = 10: 
thick solid line; z = 15: thick dashed line; z = 20: thick dot- 
dashed line), assuming a PLl spectrum for the BH emission. The 
thin solid line shows the spectrum we would obtain at z = 10, 
had we assumed that Az = in eq. (5). 



Above that the multi-color black body is exponentially cut- 
off, and the power-law component (oc E~^) emerges, as de- 
scribed in Shakura & Sunyaev 1973 (see also Salvaterra et 
al. 2005 for an application to a context similar to the one 
we are considering): 



FMC,<i{E,z) = 



Ese ^'^p 

1 __B. 

Ese 




+ A.s>E- 



0.01 eV <E <Ep 

Ep <E < lO'' eV 
otherwise. 



Ep ~ 3000 eV{MBH/MQ) 



-1/4 



(14) 



A<s, is chosen so that the energy in the power-law spectral 
component is equal to a fraction <1> of the energy in the multi- 
color black body spectral component. $ is usually taken to 
be < 1, and we wiU consider the case $ = 0.1 (MCOl SED), 
which is practically indistinguishable from all the cases with 
lower $, and quite similar to the CEise with $ = 1, too. As 
we substitute Mbh with {Mbh){z) inside eq. (14), we note 
that this spectral shape is slightly dependent on redshift^. 

For all the considered spectral shapes we chose to as- 
sume that the BH emissivity at energies below 0.01 eV or 



^ The low-energy tail of a modified blackbody spectrum is ex- 
pected to be oc E^. We neglect such slope change, as only a small 
fraction of the BH luminosity is emitted in this region of the spec- 
trum. We also note that the exponential constant was chosen to 
be 3Ep in order to ensure that the multi-color black body com- 
ponent actually has a (broad) peak at E 
Salvaterra et al. (2005). 



Early BH radiation - I: effects on the IGM 5 




Figure 2. Spectrum of the background produced by primordial 
accreting BHs, and seen by neutral-IGM baryons at z = 10. 
The four panels refer to the four BH growth history scenar- 
ios we consider (IMBH-3%: top left panel;IMBH-6%: top right 
panel; SMBH-3%: bottom left panel; BVR06; bottom right panel), 
whereas the different line types refer to three different SED for 
the BR spectra (PLl: solid line; SOSl: dashed line; MCOl: dot- 
dashed line). 

above 10® eV is 0. Such choice prevents numerical problems, 
and does not significantly affect our results. 

2.3 Results 

In Fig. 1 we show the redshift evolution of the spectrum 
of the background radiation produced by primordial BHs 
and reaching a neutral-IGM baryon. In such Figure we con- 
sider all the different BH growth histories, but only the PLl 
SED. The background level grows with time, as could be 
expected when we remember that the BH density, the av- 
erage distance among BHs, and the IGM density all evolve 
in a background-enhancing direction. In all the considered 
growth scenarios, the spectra peak at i? ~ 1 keV at all red- 
shifts: above this peak the spectrum is almost un-absorbed 
(i.e., the specific flux decreases in the same way as the in- 
put spectrum), whereas the spectrum at energies below the 
peak is shaped by the cutoff due to the IGM absorption. 
Such absorption cutoff slowly moves to lower energies. It is 
also interesting to look at the thin solid line, which illus- 
trates the effect of using Az = in eq. (5): as expected, the 
high energy part of the spectrum does not change, while the 
sharp cutoff at low energies is replaced by a much milder 
power- law decline. 

When comparing different BH growth histories in Fig. 1, 
it is clear that the normalization of the background spec- 
trum is related to the value of the product {y x pbh) at 
the relevant z. Instead, the sharpness of the low-energy cut- 
off depends on the geometrical properties of the BH spatial 



Figure 3. Redshift evolution of the total energy input per 
noutral-IGM baryon due to the background produced by primor- 
dial BHs. As explained in the caption of Fig. 2, the four panels 
refer to the four considered BH growth histories, and each line 
type refers to a different assumed BH spectrum. We also show 
the energy input from BHs with a PLl SED, had we assumed 
that Az = in eq. (5) (thin dotted line). 



distribution: in models with large values of (AIbh) the cut- 
off is very sharp, whereas it is a bit more gentle for IMBH 
models with low {Mbh}- This is important because the low 
energy part of the spectrum, albeit accounting only for a 
small fraction of the total energy in the background, is ab- 
sorbed with quite high efficiency and is a major contributor 
to the energy input e. 

In Fig. 2 we show the spectrum of the background radi- 
ation at a fixed redshift, z = 10, while varying the BH SED. 
It is clear that 'flat' (PLl, MCOl) SEDs produce larger back- 
grounds than 'steep' (SOSl) SEDs, simply because a larger 
fraction of their luminosity is emitted in the energy range 
{E ^ 100 — 1000 eV) where incomplete (or negligible) ab- 
sorption allows the build-up of the background. This is par- 
ticularly clear for the multi-component SEDs, that produce 
a quite prominent bump in a broad energy range around 
the 1 keV peak of the spectrum. Even if the peak of the 
background specific fiux is never far from ~ 1 keV, it is pos- 
sible to discern some trends: the SOSl SED tends to peak at 
slightly lower energies than the PLl SED, whereas the po- 
sition of peak of the MCOl SED depends on the chosen BH 
growth scenario, simply because its peak energy Ep depends 
on the typical BH mass. 

Finally, in Fig. 3 we show the total energy input per 
baryon as a function of redshift. Such a quantity is well 
correlated with the intensity of the background spectrum, 
especially at low energies. Thus, it increases with time, and 
the SED with the highest low energy component (MCOl) 
gives the maximum energy input. We also compare the en- 
ergy input from the reference PLl SED with that from an 
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Table 2. Fraction of the unresolved X-ray background which in our models is due to BH emission at 2 > Zdrop> in various bands. The 
numbers in parenthesis are lower limits, obtained from the assumption that the actual X-ray background is 1 — it higher than the central 
values. The cases where the emission from our models exceeds the unresolved background are in bold face. 



Model 


^drop 


u.o-z Kev 


z-o Kev 


L-Z Kev 


Z-D Kev 


0.65-1 keV<= 


IMBH-3%+PLl 


5 


0.32(0.21) 


0.19(0.10) 


0.37(0.26) 


0.42(0.13) 


n7S('n 


.065) 


IMBH-3%+SOSl 


5 


0.024(0.016) 


0.014(0.008) 


0.028(0.020) 


0.031(0.010) 


0.006(0, 


,005) 


iiViJl>Xl.-0 /O T^-lVi\_;U 1 


c: 



u.u 1 yyj.^^ J 






v.vuXjyv.yjoyj j 


0.21(0, 


,18) 


IMBH-3%-l-MC01 


6 


0.18(0.11) 


0.014(0.008) 


0.057(0.041) 


0.031(0.010) 


0.054(0, 


,045) 


IMBH-3%-l-MC01 


7 


0.046(0.033) 


0.005(0.003) 


0.016(0.011) 


0.011(0.003) 


O.Ulo(U, 


.012) 


iMBH-D7o-|-Fl-il 


5 


2.5(1.7) 


1.5(0.82) 


2.9(2.1) 


3.3(1.0) 


0.63(0, 


,52) 


IMBH-6%-fPLl 


6 


0.79(0.51) 


0.47(0.25) 


0.92(0.65) 


1.0(0.32) 


0.19(0, 


,16) 


IMBH-6%-|-PLl 


7 


0.26(0.17) 
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IMBH-6%-hSOSl 
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IMBH-6%+MC01 
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0.72(0.22) 
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0.30(0.22) 
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0.012(0, 
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BVR06-I-PL1 
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0.30(0.19) 


0.18(0.097) 


0.35(0.25) 


0.39(0.12) 
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BVR06+SOS1 
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0.014(0.007) 
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" Flux in the 0.5-2 keV band, normalized to a background level of 2.59(4.00) x 10~^ ergs~^ cm~-^ sr"^ (from the combination of the 

B04 unresolved fraction, and the Moretti et al. 2003 total X-ray background). 
Flux in the 2-8 keV band, normalized to a background level of 4.35(8.05) x 10~® erg cm~^ sr"'^, (from the combination of the B04 
unresolved fraction, and the De Luca & Molendi 2004 total X-ray background). 
Flux in the 1-2 keV band, normalized to a background level of 1.12(1.58) X 10^^ erg s^-^ cm^^ sr^-'- (from HM07). 
Flux in the 2-5 keV band, normalized to a background level of 1.31(4.26) x 10~^ ergs~^ cm~^ sr~^ (from HM07). 
« Flux in the 0.65-1 keV band, normalized to a background level of 3.28(3.94) x 10-^ ergs"! cm-2 sr"! (from HM07). 



otherwise identical model where we assumed that = in 
eq. (5). This is useful to check the eflFects of our assumption 
about A^;, and also gives us a rough estimate of the level of 
the spatial fluctuations of the energy input. The difference 
usually amounts to a factor of 2 — 3, even if it might be 
larger for the SMBH-3% and the BVR06 growth histories, 
especially at high z. 

We note that the model where the energy input from 
BHs is maximum is the one where the IMBH-6% accretion 
history is combined with the MCOl SED. lu the following, wc 
will refer to such combination as the 'extreme' model, since it 
leads to the strongest BH feedback effects (and is also close 
to the constraints from the unresolved X-ray background; 
see below). On the other hand, we will also consider the 
IMBH-3%-|-PLl model (i.e. the one combining an IMBH- 
3% history with a PLl SED) as a 'fiducial' case. 



2.3.1 Constraints from the unresolved X-ray background 

As a consistency check, wc looked at whether our models 
are compatible with measurements of the unresolved X-ray 
background from Bauer et al. (2004; hereafter B04), and 
from Hickox & Markevitch (2007; hereafter HM07). 

For such comparisou, wc obtained the spectrum of the 
background produced by the BH emission at a redshift Zdrop, 
we integrated it in the relevant energy band, and we red- 
shifted it to = assuming no absorption. 



Such a calculation implies that, at z < Zdrop, the emis- 
sivity due to BHs is 0. This is a quite crude assumption, but 
it must be remarked that observations (Steidel et al. 2002) 
suggest that the duty cycle declines with redshift (reaching 
y ~ 10"'^ at 2: = 0), and that several theoretical models in- 
clude a variation of y (e.g., in model M3 of RO04 3/ = 1 at 
z ^ 14, but y = 10~^ at z ^7). It is also possible that at 
redshifts ^ 5 — 7 an increasing fraction of the BH sources 
are detected as resolved AGN sources. We also stress that a 
fraction of IMBHs are expected to merge into larger SMBHs 
or to be ejected from the parent halos as a consequence of 
three-body encounters (see e.g. Volonteri et al. 2002, 2003). 
In these cases the IMBHs no longer contribute to the X-ray 
background. As our model does not account for these effects, 
it represents a strong upper limit for the X-ray background 
from IMBHs. 

The results of the comparison are listed in Table 2, 
where we generally adopted Zdrop — 5. It can be seen that 
two of our models (IMBH-6% growth history, combined 
with either a PLl or a MCOl SED) exceed the observed 
background in at least one band. For such cases (and also 
for the IMBH-3%-|-MC01 case, where the contribution from 
BHs exceeds half of the unresolved X-ray background in the 
0.5 — 2keV band) wc also list the result wc would obtain 
with Zdrop = 6 or 7, which clearly show that the constraints 
from the X-ray background can be easily satisfied also by 
these models, provided that Zdrop ~ 6. Thus, we note that 
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Figure 4. Redshift evolution of the hydrogen ionization fraction 
Xff. The order of the panels and the meaning of the various line 
types are the same as in Fig. 2, except for the thin dotted line, 
which represents the ionization evolution in a model without any 
BH emission. 



Figure 5. Redshift evolution of the IGM temperature in re- 
gions outside ionized bubbles. The order of the panels and the 
meaning of the various line types are the same as in Fig. 2, ex- 
cept for the thin dotted line, which represents the temperature 
evolution in a model without any BH emission, and for the thick 
dotted line representing the CMB temperature. 



the choice of Zdrop (and, in general, the fate of IMBHs in 
the lower redshift range we consider) is quite crucial for our 
models. In the rest of this paper we will use Zdrop = 5 for 
all the models, and our plots will extend to such redshift. 



3 INFLUENCE ON THE IGM EVOLUTION 

We looked at the effects of the energy input due to the back- 
ground radiation produced by primordial BHs on the ther- 
mal and chemical evolution of the IGM. We employed a sim- 
plified version of the code described in Ripamonti, Mapelli 
& Ferrara (2007; hereafter RMF07; but see Ripamonti et al. 
2002, and Ripamonti 2007 - hereafter R07 - for more detailed 
description of this code), in order to look at the evolution 
of the IGM under the infiuence of the energy input we cal- 
culated in the previous Section. 

Such a code follows the gas thermal and chemical evo- 
lution. The chemistry part deals with 12 chemical species 
(H", H+, H", D°, D+, He°, He+, He++, Ha, H+, HD, and 
e~ ) , and includes all of the reactions involving these species 
which are listed in the Galli & Palla (1998) minimal model 
for the primordial gas, plus some important extension (e.g., 
it considers the ionizations and the dissociations due to the 
energy input we are introducing) . The thermal part includes 
the cooling (or heating, if the matter temperature is lower 
than the CMB temperature) due to molecules (H2 and HD), 
to the emission from H and He atoms, to the scattering of 
CMB photons off free electrons, and to bremsstrahlung radi- 
ation. Furthermore, it accounts for the cooling/heating due 



to chemical reactions, and for the heating due to the energy 
input we are considering®. 



3.1 IGM ionization and temperature 

Figs. 4 and 5 show the effects of the BH emission upon the 
ionization level (in particular, the hydrogen ionized fraction 
xh = n(//+)/[n(ff°) -I- n{H'^)] ) and the temperature of 
the IGM Tfe (all these quantities are calculated outside the 
ionized bubbles close to radiation sources) . In all the models 
we consider, the BH emission starts altering the neutral IGM 
at 2 ~ 15 — 20. After that there is a steady increase in both 
Xh and Tk- The increase of xh stops only when the IGM is 
completely ionized (however, such condition is reached only 
in the most extreme of our models, and only at a redshift 
~ 6). Instead, Tk stops increasing once it reaches a level 
(~ lO** K), where atomic cooling is important: in the models 
where BH emission is assumed to be strongest (IMBH-6% 
with MCOl spectrum) this happens at 2; ~ 10, but 2 6 — 8 
is a more typical range. 

It must be noted that in the lower redshift range we 
consider (say z ^ 10) our models start suffering from several 



® Other than introducing the energy input as calculated in the 
previous Section, the code differs from the version described in 
RMF07 because we introduced the cooling through He lines and 
bremsstrahlung (the rates were taken from Anninos et al. 1997), 
and we splitted the energy input into the heating, ionization and 
excitation components by using the expressions given in eq. (11), 
rather than the Chen & Kamionkowski (2004) approximations. 
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problems. First of all, the energy input we employ is calcu- 
lated for a neutral medium, whereas in some of our models 
xh ^ 0.5 already at 2 ~ 7 — 8. Then, we are overestimat- 
ing the energy input ^. Second, we are assuming that the 
IGM density remains constant at its average unperturbed 
cosmological value, whereas this approximation becomes in- 
creasingly problematic as structures start to form. Third, 
there might be some level of metal enrichment (altering both 
the heating and the cooling rates) even in regions which 
are far away from the most luminous sources. Finally, we 
are completely neglecting the contribution to heating and 
reionization which is due to stars, which is likely to be sub- 
stantial at relatively low redshifts. However, our calculations 
should still be reasonably accurate until the end of the so- 
called 'overlap' phase of reionization (probably not far from 
2^7 — 8, see Section 5.1), provided that they are taken to 
represent conditions in regions which were not yet ionized. 



3.2 CMB angular spectrum 

The cosmic heating and the contribution to reionization due 
to BHs might also leave some imprint on the CMB spectra. 
In order to study this effect, we implement ionization and 
gas temperature evolution due to BHs in the version 4.5.1 
of the public code CMBFAST (Seljak & Zaldarriaga 1996; 
Seljak et al. 2003). 

Fig. 6 shows the temperature — temperature (TT) , po- 
larization — polarization (EE) and temperature — polar- 
ization (TE) spectra of the CMB in the case in which the 
contribution from BHs is accounted for. In particular the 'ex- 
treme' case (IMBH-6%+MC01, solid line) and the 'fiducial' 
one (IMBH-3%-|-PLl, dashed line) are shown. They are also 
compared with the spectra obtained without contributions 
from stars and/or BHs (thin dot-dashed line) and with the 
spectra derived assuming Thomson optical depth Te ~ 0.09 
and a sudden reionization at 2; ~ 11, consistent with the 3-yr 
Wilkinson Microwave Amsotropy Probe ( WMAP) data. 

As one can expect, no significant differences appear be- 
tween the four cases in the TT spectrum. The 'fiducial' and 
the 'extreme' BH model differ from the thin dot-dashed line 
both in the EE and in the TE spectra, at low multipoles 
{I ^ 20). However, the contribution of BHs to the TE and 
EE spectra, even in the extreme case (IMBH-6%-|-MC01), 
is smaller than (or comparable to) the best fit of the 3-yr 
WMAP data (dotted line). Thus, all the scenarios consid- 



The situation is actually quite complicated. The reduction in 
the heating rate is slower than what could naively be expected 
[e oc (1 — xh)~^] from the increase of xh, because the bulk 
of the cross section is due to He, which is harder to strip of 
its electrons (see e.g. Thomas & Zaroubi 2007). For example, 
in the 'extreme' model, at redshift 6 Xf{ ~ 0.9, but a;^e+-|- = 
n{He++)/[n{He°) + n(He+) + n{He++)] ~ 0.1, and by using 
eq. (2) we are overestimating the energy input e only by a fac- 
tor of ~ 2.5 rather than (1 — xh)~^ ~ 10. Furthermore, the 
assumption that the IGM is completely neutral also leads to an 
overestimation of the optical depth r, and the background radia- 
tion (and energy input) is correspondingly underestimated. Then, 
our energy input rates are essentially correct, except for our 'ex- 
treme' model at 2; ^ 7, where we might be overestimating e by a 
factor of 2. 
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Figure 6. Effects of BH emission upon the CMB angular spectra. 
Temperature-temperature (top panel), polarization-polarization 
(central panel) and temperature-polarization (bottom panel) 
spectra are shown. Thick dotted line: CMB spectra derived as- 
suming Thomson optical depth Te = 0.09 and a sudden reion- 
ization model (consistent with the 3-yr WMAP data,); thick solid 
and dashed lines; CMB spectra derived assuming energy injection 
from the BHs in the 'extreme' (lMBH-6%+MC01, i.e. the case 
where BH energy input is strongest) and the 'fiducial' (IMBH- 
3%-|-PLl) cases, respectively; thin dash-dotted line: CMB spectra 
derived assuming no reionization and no contribution from BHs. 



ered in this paper (even IMBH-6%+MC01) do not violate 
the limits posed by CMB observations. 

Furthermore, the Thomson optical depth which can be 
directly derived from the ionization history shown in Fig. 4 
is Te < 0.07 (re = 0.027 and 0.064 in the 'fiducial' and 'ex- 
treme' case, respectively), smaller than the best fit to 3-yr 
WMAP data (re = 0.09 ± 0.03, Spergel et al. 2007). Thus, 
our BHs might give a partial contribution to the reioniza- 
tion, but are not its exclusive source, in agreement with 
previous work (e.g. RO04; Ricotti et al. 2005; Z07+). 



3.3 21 cm radiation 

3.3.1 Basic definitions 

The spin temperature of the 21-cm transition can be written 
as (see e.g. Field 1958, 1959; Kuhlen, Madau & Montgomery 
2006; Valdes et al. 2007; Z07+) 

^ _ r» + TCMB + (Vk + ya)Tk , , 

J- spin — 111 ' \'-^) 

where T, = 0.068 K corresponds to the 21-cm transition en- 
ergy, TcMB is the CMB temperature, is the IGM kinetic 
temperature, and y^ and ya are the kinetic and Lyman a 
coupling terms, respectively. 

The kinetic coupling term is 
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Figure 7. Redshift evolution of the neutral H spin temperature. 
The order of the panels is the same as in Fig. 2; the dotted line 
represents the CMB temperature, thick lines represent the spin 
temperature, and thin lines represent the IGM temperature. Con- 
tinuous lines refer to models with a PLl SED, dashed lines to 
models with a SOSl SED, and dot-dashed line to models with a 
MCOl SED. 



Vk 



(CH+Ce + Cp), 



Alo Ik 

where Aw ~ 2.85 x lO"'^^ s^^ is the Einstein spontaneous 
emission rate coefiicient (Wild 1952), and Ch, Ce and Cp 
are the de-excitation rates due to neutral H, electrons and 
protons, respectively. They are given by the fitting formulae 
from Kuhlen et al. 2006 (see also Field 1958, 1959; Smith 
1966; Allison & Dalgarno 1969; Zygelman 2005): 



Ch 



3.1 X 10 



He 7e 

3.2 {np/nu) Ch, 



(18) 

(19) 
(20) 



where uh, Ue and Up are the neutral H, electron and proton 
number densities, and 



logi 



1 cm3 s~i 



-9.607-H0.51ogio exp 



1 / T 
1800 O^Sioj^ 



The Lyman a coupling term is given by 



167r r, vre^ 



/i2 Jo, 



(22) 



27Aio Tfc me c ' 

where e and me are the electron charge and mass, /12 — 
0.416 is the oscillator strength of the Lyman a transition, 
and Jo is the intensity of Lyman a photons which are due 
to coUisional excitations from thermal electrons, to hydro- 
gen recombinations, and to coUisional excitations from X-ray 
energy absorption. Jo is then 



Jo 



he 



4ttH{z) 



rlerlH'ye.H + UeUpa'^lp + 



hUc, 



(23) 



Figure 8. Redshift evolution of the brightness temperature dif- 
ference with respect to the CMB <5T(,. The order of the panels 
and the meaning of the various line types are the same as in Fig. 
2, except for the dotted line, which represents 5T), for a model 
without any BH emission. 



where ^e,H =i 2.2 x iQ-Sg-(iiS400K)/Tk cm^s'^ is the col- 
lisional excitation rate of neutral H atoms by electron im- 
pacts, Va ~ 2.46 X 10^^ Hz is the Lyman a frequency, and 
C('^2p is the effective recombination coefficient to the 2^P 
level (including recombinations to the 2^P level, plus re- 
combinations to higher levels that end up in the 2^P level 
through all possible cascade paths). We adopted a simple 
fit to the Pengelly (1964) results for ctj^p' assuming case A 
recombinations* : 



iff 

^22p 



{Tk) 



1.67X lo-"^^°■^l-(^/"^)'°^^^^ 



(24) 



where T4 =Tfe/(10^ K). 

Once the spin temperature is known (from eq. 16), it is 
convenient to express the resulting 21-cm radiation intensity 
as the differential brightness temperature between neutral 
hydrogen and the CMB, which is an observable quantity: 



erri ^ Tspin TCMB , c \ 

(1 + ^) "^^(1+^-)' 



(25) 



where 5p = (p — p)/p is the cosmological density contrast 
in the considered region (here we will consider only the case 
5p = Qi), and T21 is the IGM optical depth at an observed 
wavelength of 21(1 -f z) cm, 



T21 



riH 



32-K kBl'2lH{z) Tspin 



(26) 



where ks is the Boltzmann constant, and 1^21 — 1.421 x 
10^ Hz is the (rest-frame) frequency of the 21-cm line. 



Results for case B recombinations differ only slightly. 
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3.3.2 Results for pure BH coupling 

Fig. 7 shows the redshift evolution of Tspin, under the as- 
sumption that only the radiation produced by BHs is im- 
portant. 

In all these models, Tgpin remains very close to Tcmb 
(and to the predictions of models with no BH emission) un- 
til a ^ 9 — 15, i.e. until Tk finally becomes much larger than 
Tcmb- After that, the difference between Tspin and Tcmb 
becomes significant, and in models with strong BH emission 
it can amount to 90 K. Apart from the amplitude of this 
maximum difference, the strength of BH emission also influ- 
ences the redshift when it is reached: in models with weak 
BH emission (e.g. most models where a SOSl SED is as- 
sumed), Tspin— Tcmb keeps increasing, and is largest at the 
lowest considered redshift (though this maximum is quite 
low); whereas in models with strong BH emission (e.g. all 
those where a MCOl SED is assumed) Tspi„ — TcMB reaches 
a relatively high maximum at z ~ 6 — 9, slowly decreasing 
afterwards. The main reason is that in models with high BH 
emission the ionized fraction easily reaches the regime (at 
xh ^ 0.1) where fexc (and Jo and ya with it, as Jo is domi- 
nated by the term due to coUisional excitations from X-ray 
absorption) starts dropping very fast, rather than being ap- 
proximately constant (see fig. 4 of Shull & Van Steenberg 
1985). 

In Fig. 8 we show the corresponding evolution of the 
differential brightness temperature STi,. Such evolution es- 
sentially mirrors the one of Tspin — Tcmb : it remains close to 
until z ~ 9 — 15, and then starts growing, reaching maxima 
between ^ 5 and ^ 18 mK, depending on the strength of 
the BH emission. Again, in models with weak BH emission 
the maximum is reached at the lowest considered redshift, 
whereas in the other models it is reached at 2 ~ 6 — 9. The 
main difference with the evolution of Tspin — Tcmb is that 
the decline after the maximum is faster, since the high IGM 
ionization level in models with strong BH emission reduces 

also T21. 

It must be remarked that such an evolution of 5Tb in 
the neutral patches of the Universe at z ^ 12 should be de- 
tectable with the new generation of radio experiments, such 
as LOFAR, MWA, 21CMA and SKA^. For example, LOFAR 
will probe the 21 cm emitted from the IGM in the redshift 
range of 6-11.5 and will be sensitive to scales from a few 
arcminutes up to few degrees and will be able to statisti- 
cally detect the 21 cm brightness temperature down to «5 
mK (de Bruyn, Zaroubi & Koopmans 2007, de Bruyn et al. 
2008). However, we stress that these effects are observable 
only before the end of the reionization epoch (see Section 
5.1). 



for example, most models of reionization (e.g. RO04, and 
references therein) assume that the stellar contribution was 
dominant over the one from BHs. This is supported by ob- 
servations of the unresolved X-ray background, whose level 
is difficult to reconcile with the hypothesis that reionization 
is due to BH emission (e.g. Dijkstra et al. 2004). Further- 
more, our models do require the presence of stellar radiation, 
as even the 'extreme' one (IMBH-6%-|-MC01) is unable to 
reionize the Universe before « ~ 6, and is therefore incom- 
patible with observations of quasars and Lyman a emitters 
at « - 6 - 7 (Becker et al. 2001; Djorgovski et al. 2001; Fan 
et al. 2001, 2002, 2003; White et al. 2003; Kashikawa et al 
2006; lye et al. 2006; Ota et al. 2007). 

Since we are looking at the evolution of the IGM in re- 
gions which are quite removed from BHs (and, consequently, 
from the bulk of stellar emission) and are reionized late'^'^, 
the omission of the stellar contribution from our calculations 
is mostly justified. In fact, the neutral IGM we are consid- 
ering is almost perfectly transparent to radiation with fre- 
quencies below the H ionization threshold (13.6 eV): such 
photons can travel cosmological distances, but are unable to 
significantly affect the IGM. On the other hand, the ioniz- 
ing photons emitted from stars arc typically absorbed at the 
edge of the ionized regions around stellar sources, since they 
are not energetic enough to cross significant distances in a 
neutral IGM: for this reason, their effects are purely local. 

Lyman a photons are the only relevant exception. In 
fact, although the Lyman a cross section is very high, such 
photons can scatter many times before exiting the reso- 
nance; more importantly, the redshifting of photons with 
energies slightly higher than 10.2 eV 'into' the resonance 
ensures a roughly uniform Lyman a radiation field also in 
neutral regions. 

Ciardi & Salvaterra (2007; hereafter CS07) found that 
the Lyman a radiation field can moderately heat the IGM: 
the heating rate taken from their models dominates over 
that of all our models at z ^ 15, and takes the IGM tem- 
perature to ^ 30 K at z ~ 15 — 20. On the other hand, at 
« ^ 10 the Lyman a heating rate should be much smaller 
than those of our models^^. More importantly, CS07 find 
that for z ^ 27 the intensity Ja,, of the Lyman a back- 
ground is much higher than the level [Ja, coupling ^ 10^ (l-(- 
z) erg cm^"^ s^^ Hz^^ sr~^] which should couple Tspin to Tj, 
rather than to Tcmb (see Ciardi & Madau 2003). 

In our case, it is reasonable to neglect the heating effects 
of the Lyman a background, although this will lead us to 
somewhat underestimate the IGM temperature at z ^ 15. 
But it is very important to add the effects of the Lyman a 
background to our estimation of Tspin (and 5Tb). 

This can be done very easily by modifying eq. (16) into 



3.3.3 Results for BH and stellar coupling 

In the previous Subsection we considered the evolution of 
the 21-cm emission under the effects of BH emission only. 

But it is largely believed that stellar emission played a fun- 
damental role in the evolution of the primordial Universe: 

® http://www.lofar.org 

http : // www . haystack, mit . edu/ ast /arrays /mwa / 
http: / /21cma. bao.ac.cn/ 
http: / /www.skatelescope.org 



It is natural to wonder down to which redshift such neutral re- 
gions actually exist. Here, we wiU simply assume that they survive 
down to 2 ~ 5, and look into their properties. Such hypothesis 
will be discussed in Section 5.1. 

The plots in CS07 actually stop at z = 10; but it is pretty 
clear that in their model the IGM temperature is growing at a 
much slower rate than in our models. It is also worth noting that 
some of the CS07 assumptions (e.g. the values of the parameters 
fgas and /*) are quite extreme, and would result in a very early 
complete reionization. More realistic assumptions would result in 
a significant delay in the rise of Tgpin ■ 
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Figure 9. Redshift evolution of the neutral H spin temperature, 
when the Lyct coupling due to stellar radiation (but not the stellar 
radiation heating effects) is kept into account. The order of the 
panels is the same as in Fig. 2. The dotted line represents the 
CMB temperature, thick lines represent the spin temperature, 
and thin lines represent the IGM temperature. Continuous lines 
refer to models with a PLl SED, dashed lines to models with a 
SOSl SED, and dot-dashed Une to models with a MCOl SED. 



Tspin — 



r* + TcMB + iVk +ya + ya,*)Tk 

1 + yk + Va + ya,. 



(27) 



where ya,, accounts for the additional coupling due to the 
Lyman a background due to the stars^'^, and is approxi- 
mately given by (see CS07): 



10^ 



AioTk 

After approximating Jq,, with the expression 



(28) 



Ja,*{z) 



erg/ (cm^ sHz sr) 




2 > 30 

30 > z > 10 (29) 
z < 10, 



which is a moderate underestimate of the Ja,* curves shown 
in fig. 1 of CS07, we have recalculated the evolution of Tspin, 
and 5Tb. The results are shown in Figs. 9 and 10. 

In this case, Tspin (Fig. 9) is almost perfectly coupled to 
the kinetic temperature, and the difference Tspin — Tcmb — 
Tspin easily reaches the 10^ - lO"* K range. Also 5Tb (Fig. 10) 
is affected. Here we focus on the relatively low redshifts 
which will be explored by 21-cm experiments (e.g. LOFAR), 
where the effects of BH emission lead to differential bright- 
ness temperatures which can reach 20 — 30 mK at redshifts 
~ 8 — 15. Instead, at high redshifts (say, z ^ 15) 5Tb can 
reach very high negative values (in the —200 to —300 mK 

^■^ Also BHs produce a Lyman a background; but its intensity 
is much lower than the one due to stars, and the corresponding 
coupling term is always much smaller than ya- 



Figure 10. Redshift evolution of the brightness temperature dif- 
ference with respect to the CMB 5Tb, when the Lya coupling due 
to stellar radiation is kept into account. The order of the panels 
and the meaning of the various line types are the same as in Fig. 
2. The insert in the bottom-right panel shows the same quanti- 
ties (for the BVR06 case; the other cases are qualitatively similar) 
on a much wider 5Tb scale. The dotted line stopping at 2 = 10 
comes from Fig. 5 (bottom panel, solid line) of CS07, and shows 
5Tb for a model with stellar Lyman a coupling and heating, but 
no BH emission. This line represents an upper limit on 5Tb in the 
absence of BH heating. 



range); but in such redshift range the results of CS07, pre- 
dicting a minimum value of STb ~ —170 mK at z ^ 24 are 
likely more correct because they include also the the heating 
effects of the stellar Lyman a background. 

We point out that our results, especially those about 
5Tb, depend only weakly from the very high level of Ja,* 
given in the CS07 paper: the effects of lowering Ja,* to 
a more realistic level, e.g. a fraction 0.1 (or even 0.01) of 
the amount given by eq. (29) are a certain reduction (from 
~ 8-10000 K to 3-5000 K) of the level where Tsp^n reaches 
a low-redshift 'plateau', and a much smaller change in the 
evolution of 5Tb. Then, our predictions about 5Tb observa- 
tions are quite independent from the assumptions of CS07. 
Instead, for the model where no BH feedback is included, a 
reduction by a similar factor in the Lyman-a heating rate 
in the CS07 models would result in a much lower 5Tb value 
than shown by the dotted curve in Fig. 10. 



4 INFLUENCE ON STRUCTURE 
FORMATION 

In the previous section we have shown that the energy input 
from BHs can substantially heat the IGM. In turn, this is 
likely to affect the formation of galaxies: as the cosmologi- 
cal Jeans mass depends on T^''^, the baryonic component of 
small fluctuations might become unable to collapse and form 
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stars because of the temperature increase. But the effects of 
BH radiation are not limited to the heating, since the in- 
crease in the H ionized fraction also enhances the formation 
of H2, which is the most important coolant in metal- free 
gas at temperatures 10'* K: such enhancement would fa- 
cilitate the formation of stars within small halos. Then, we 
investigated the influence of BH energy input on structure 
formation with a method which accounts for such competing 
effects, and which was already employed in the RMF07 and 
R07 papers. 

We used the full code (instead of the simplified version 
used in the previous Section) described in Section 3, in or- 
der to follow the evolution of spherically symmetric halos 
of different masses, virializing at different redshifts. Such 
evolution took into account all the physics included in the 
simplified version we already described, plus the treatment 
of gravity, of the hydrodynamical evolution of the gas, and 
of the dissociation of H2 molecules due to Lyman- Werner 
(11.2 eV <hv < 13.5 eV) photons emitted by BHs". Dark 
matter (DM) gravitational effects are included as described 
in Sec. 2.1.3 of R07: the DM final density profile is assumed 
to be a truncated isothermal sphere with ^ = 0.1 (i.e., the 
final core radius is assumed to be 1/10 of the virial radius). 



4.1 Critical mass 

As in the RMF07 and R07, we classified halos as collaps- 
ing if they reach a maximum density larger than pcou = 
1.67 X 10~*® gcm~^ ~ lO^mff cm~^ (a value high enough to 
suggest that the formation of a luminous object is well un- 
der way) in less than an Hubble time after their virialization 



(at Zvir), i.e. at a redshift z^. 



[0.63(1 Z„ir)] - 1. 



This classification criterion is roughly comparable to the 
collapse criterion of Tegmark et al. (1997): in analogy with 
such paper (and with RMF07 and R07), we define the crit- 
ical mass Mcrit[zvir) as the minimum mass of a collapsing 
halo virializing at Zvir- 

In Fig. 11 we compare the evolution of Merit which is 
obtained for each of our BH models with the same evolution 
in the unperturbed (e = at all redshifts) case, and with 
the evolution of the mass 



MH{Zv^r) ^ 1.05 X 10" Mq (1 + Z^.r)'^^^ 



(30) 



of halos with a virial temperature Tvir = 10* K (assuming 
a mean molecular weight /i = 1.23, as appropriate for a 
neutral medium), above which the cooling due to atomic H 
becomes dominant. 

The BH energy injection has negligible effects upon 
Merit for Zvir ~ 15, but its effects become increasingly im- 
portant at later times: at « = 10 the BH energy input in- 
creases Merit by a factor between 1.8 (SMBH-3%+S0Sl 
model) and 40 (IMBH-6%-FMC01 model). At lower red- 
shifts the BH effects are even larger: in the models with the 
strongest BH feedback. Merit can become Si Mr, although 



For this last effect we use the reaction rate given by Abel et al. 
1997 (reaction 27); the flux of photons at 12.87 eV was obtained 
through the formalism described in Section 2.1, but assuming that 
photons at frequencies corresponding to the lines of the Lyman 
series of hydrogen were completely absorbed. No stellar emission 
was assumed. 




20 15 12 10 8 20 15 12 10 
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Figure 11. Redshift evolution of the critical mass. The order of 
the panels and the meaning of the various line types are the same 
as in Fig. 2, except for the thin continuous line, which represents 
^erit for a model without any BH emission, and for the dotted 
line, representing the mass Mh of halos with virial temperature 

Tv^r = 10* K. 



the onset of atomic cooling slows down the increase of Merit- 
in such models, BH feedback prevents the formation of stars 
inside mini-halos cooled by molecules at z 9. 



4.2 Gas retention 

RMF07 suggested that one possible feedback effect of the 
energy input from decaying/annihilating DM particles is to 
reduce the amount of gas which actually ends up within the 
potential wells of virialized halos. As the feedback effects 
of BHs are much stronger than those of DM decays and 
annihilations, we looked at whether these same effects are 
important in our simulations. 

To this purpose, we define fret as the ratio of the mass 
of gas which is retained inside the virial radius of a halo (at 
the time when our simulations are stopped) with respect to 
the baryonic mass expected from cosmology. For a halo with 
total mass M^aio and baryonic mass Algas'- 



Mgas n 

n 



Mhalo 



(31) 



In Fig. 12 we show the dependence of fret upon Mhaio, 
for halos virializing at z^ir ~ 10 and for all our BH models, 
plus the unperturbed case. 

Generally, models with BH feedback exhibit a sharp 
transition at Mhaio — Merit, going from fret ~ to 
fret ~ 0.7, whereas in the unperturbed case fret increases 
quite smoothly from ~ 0.4 to ~ 0.9 (this last value is in good 
agreement with numerical simulations by Grain et al. 2007). 
This threshold effect is due to the hydrodynamical effects of 
the BH heating, combined with the depth of the DM poten- 
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Figure 12. Fraction of gas retained by a halo at the end of our 
simulations in halos virializing at z^ir = 10, as a function of the 
total halo mass. The order of the panels and the meaning of the 
continuous, dashed and dot-dashed lines are the same as in Fig. 2; 
the dotted lines refer to the unperturbed model. The thickness of 
the line indicates whether a certain halo mass is below (thin line) 
or above (thick line) the critical mass. The seesaw behavior for 
high masses is purely numerical (i.e. due to the discrete number 
of shells). 



Figure 13. Central gas temperature at the end of our simulations 
in halos virializing at z„ir = 10, as a function of the total halo 
mass. The order of the panels and the meaning of the continuous, 
dashed and dot-dashed lines are the same as in Fig. 2; the dotted 
lines refer to the unperturbed model. The thickness of the line 
indicates whether a certain halo mass is below (thin line) or above 
(thick line) the critical mass. 



tial well. In fact, the heating induced by the BHs amplifies 
the pressure gradients and tends to prevent the gas from 
falling inside the DM potential well. If the DM potential 
well is below a certain critical value (i.e. if Mhaio ^ Merit), 
the final gas over-density is generally ^ 10, whereas the DM 
over-density is ^ 1000. On the other hand, if the gravity of 
a halo is strong enough (i.e. if Mhaio ~ Merit), the heat- 
ing induced by the BHs cannot counteract the gravitational 
pull, and the halo will retain most of its gas, which will cool, 
collapse and form luminous objects. 

This is confirmed by Fig. 13, where we show the final 
temperature of the gas at the centre of halos virializing at 
Zvir = 10, as a function of the halo mass. In halos with 
mass < Merit the gas temperature is much higher, if BH 
heating is present, than in the unperturbed case, and it is 
close to the temperature Tk of the IGM (1000 - 10000 K). 
For masses > Merit the final temperature in presence of BH 
heating is similar to the unperturbed case (200 — 400 K, 
much lower than the temperature of the surrounding IGM), 
as the gas in the centre of the halo was able to condense 
and cool. The transition in the case of the gas temperature 
is even sharper than in the case of fret, probably because 



the density dependence of the cooling rate will lead to a 
'runaway' cooling as soon as the density starts to increase. 

Fig. 12 also shows that in models with strong 
BH feedback (IMBH6%-hPLl, IMBH6%-f MCOl, and 
BVR06-I-MC01) the transition from fret ~ to fret ~ 0.8 
is not as sharp as in the other cases we consider (remark- 
ably, such difference is not present in the temperature plots 
of Fig. 13). In such models, halos with masses in the range 
Merit < Mhaio ~ SMerit are relatively poor in gas, despite 
being able to form luminous objects at their centre. Such a 
luminous but gas-poor halo population starts developing at 
z ~ 12, and becomes increasingly important when lower red- 
shifts are considered: for instance, at 2 ~ 8 this population is 
present also in models with intermediate BH feedback, and 
can span a factor of ~ 10 in mass. If such objects actually ex- 
ist and survive until present, they should be characterized by 
a high M/L ratio, a low gas content, and a mass ~ 10® Mq. 
Such properties remind us of the dwarf spheroidal galaxies 
of the Local Group (see Mateo 1998), although it might just 
be a coincidence. Further investigation is needed to address 
this issue. 



In most of the regimes we are considering, the cooling rate is 
due to H2 molecules. An increase in density results in both an 
increase in the cooling rate per molecule (which is oc p) and an 
increase in the abundance of molecules. 
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5 DISCUSSION 

5.1 Do neutral regions exist below redshift 11? 

The 3-yr WMAP results (Spergel et al. 2007) for the electron 
scattering optical depth can be interpreted as indicating a 
sudden reionization at z ~ 11. In such a scenario, neutral 
regions essentially cease to exist as soon as reionization hap- 
pens, and the effects of BHs at ^ 11 would become negli- 
gible: 

• The 21 cm brightness temperature differences would be 
severely quenched because of the lack of neutral H {5Tb 
is proportional to the density of neutral H atoms). Even if 
it were not, it is reasonable to expect that the reionization 
takes the IGM temperature Tk to ~ lO'' K, and the BH heat- 
ing would be unable to drastically change T^; the changes 
in Tspin would be even smaller. 

• The high IGM temperature we just mentioned would 
probably have important effects on structure formation; but 
that is a feedback effect from stellar sources, rather than 
from BHs. 

In short, the effects of BHs can be clearly observed only 
at redshifts before the end of reionization process. In the 
case of a sudden reionization at 2 ~ 11 they would become 
extremely difficult to detect, except perhaps in our models 
with the strongest BH feedback. 

However, the sudden reionization scenario appears un- 
realistic. In fact, practically all the theoretical models pre- 
dict that the reionization process is quite extended in time. 
In particular, the most recent numerical simulations (e.g. 
Iliev et al. 2007; Santos et al. 2007; Zahn et al. 2007; 
Mesinger & Furlanetto 2007) essentially agree in the pre- 
diction that the end of the overlap phase (i.e., the time 
when the volume filling factor of neutral regions becomes 
negligible) is at redshift 6.5 ^ Zoveriap ^ 8. Fig. 3 of Santos 
et al. 2007 is particularly useful for our purposes, since it 
includes not only the evolution of the volume-averaged ion- 
ization fraction (solid line), but also a similar curve where 
complete ionization is assumed within the ionized 'bubbles' 
(dashed line); it is quite reasonable to expect the volume 
filling- factor of neutral regions to drop below ^0.1 when 
such curve exceeds 0.8 — 0.9, i.e at z ^ 7 — 8. 

Such behavior is broadly consistent also with analytical 
models such as the one presented in Choudhury & Ferrara 
2006 (in their fig. la the volumc-avcragcd neutral fraction 
goes below 0.1 already at « ~ 9, but declines below 0.01 
only at z ^ 6). It is also important to point out that such 
low values for Zoveriap arc usually obtained in models based 
on the 3-yr WMAP data, whereas models based upon 1-yr 
WMAP data (Kogut et al. 2003; Spergel ct al. 2003) lead to 
significantly higher values of Zoveriap (for instance, see Iliev 
et al. 2007, which presents the results of simulations based 
on both sets of parameters). 

Furthermore, observations of the Gunn-Peterson 
troughs in the spectra of quasars at z ^ 6 (sec e.g. Fan et. al. 
2002; Fan et al. 2006), and measurements of the evolution 
of the density of bright Lyman a emitters at 5.7 ^ z ^ 7 
(Kashikawa ct al. 2006; lyc et al. 2006; Ota ct al. 2007) 
might hint that we are actually observing the final stages of 
overlap; but the interpretation of the data is difficult and 
the issue is still under intense debate (see e.g. Malhotra 



& Rhoads 2004; Fan, Carilh & Keating 2006; Mesinger & 
Haiman 2007; Dijkstra, Wyithe & Haiman 2007). 

In short, it seems reasonably likely that the volume fill- 
ing factor of neutral regions remained significant ( ^1 0.1) at 
least until redshift 7 — 8, and maybe even at lower redshifts; 
nonetheless, it is also possible (e.g. if the 3-yr WMAP results 
are underestimating Te) that a dearth of neutral regions at 
z ^ 10 will prevent the detection of the effects we discuss. 

5.2 Comparison with previous works about BH 
feedback 

Our analysis of feedback effects from high-rodsliift BHs has 
many links with the former study by RO04 (and also with 
Ricotti et al. 2005). However, there are some crucial differ- 
ences between the assumptions in the two sets of models, 
which lead to important differences in the results. 

• The most important difference is likely to be in the 

growth histories of the BH densities. All of the RO04 mod- 
els reach psH ^ 10^ M© Mpc^"^ (in their notation, lubh ~ 
1.7 X 10~^; cfr. the lower panel of their fig. 2) at red- 
shifts > 15, whereEis at z = 15 none of our models exceeds 
Pbh = 10^ Mq Mpc~^. This difference becomes less impor- 
tant when going to lower redshifts. The IMBH-6% growth 
history actually overtakes the RO04 predictions at « ^ 6 — 7. 
But all the other growth histories we consider are at most 
comparable to the RO04 models even at 2 = 5. 

• We assume a constant duty 
cycle (y = 0.03, 0.06, or 0.10, depending on the model), 
whereas in the RO04 models this quantity strongly depends 
on redshift (see the bottom panel of their fig. 3): it is as- 
sumed to be 1 at high redshifts {z > 14, 2 > 19, or 2 > 24), 
but rapidly falls to 10~^ when lower redshifts are considered 
(z < 13, OT z < 8). 

• RO04 restrict their analysis to an intrinsically absorbed 
Sazonov et al. (2004) spectrum; their treatment of radiation 
transfer is more detailed than in the present paper, but as 
they are not limiting themselves to the neutral-IGM, their 
background spectrum is likely to extend to lower energies 
than ours, resembling the thin solid line in Fig. 1. 

• The RO04 models include also a stellar contribution. 

Because of all these differences, the RO04 models pre- 
dict a much larger energy injection into the IGM (at z ^ 15 
the difference can easily amount to a factor of ^ 10^). At 
lower redshifts (2^8 — 9) such difference is erased (or even 
reverted), mostly because of the reduction of the RO04 duty 
cycle. 

Taking into account these differences, the results of the 
current paper are reasonably consistent with those of RO04. 
In fig. 5 of RO04, the ionized fraction and the IGM temper- 
ature are shown for different models. Complete ionization 
is achieved already at 2 ~ 7 — 8, while in our models xui 
is always less than 1 at 2 > 6 (but in most of our models 
complete reionization is never reached). This difference is 
simply explained by the presence of a stellar component in 
the RO04 semi-analytical model. For the same reason, IGM 
temperatures of 10* K are reached at 2 ~ 8 — 10 in our paper 
and at 2 20 - 25 in RO04. The Thomson optical depth 
derived by RO04 is 0.1 < Te ^ 0.2, but a fraction Te ~ 0.06 
is due to stars. Thus, the contribution of BHs to the Thom- 
son optical depth in RO04 models is Te « 0.04 — 0.14, 
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Figure 14. Redshift evolution of tlie total energy input per 
baryon due to the background produced by X-ray emission as- 
sociated with star formation (thick lines; the two lines refer to 
power-law SEDs with different photon index: F = 1.5 for the 
dashed line, F = 2.0 for the solid line, and F = 2.5 for the dotted 
line), assuming fx = 1. For comparison, we show also the energy 
input in models where BH emission was considered (IMBH-6% 
-I- MCOl: thin dot-dashed Une; IMBH-3%-|-PLl: thin solid hue; 
BVR06 + SOSl: thin dashed line). 

which is consistent with our findings (te 0.07). Further- 
more, RO04 aim to reproduce the Thomson optical depth 
derived from 1-yrWMAP results (re ~ 0.17±0.05), which is 
considerably higher than in the 3-yr WMAP measurements 
(re ~ 0.09 ±0.03). 

Ricotti et al. (2005) also study the effects on the 21- 
cm line; but, in their fig. 10, (5Tb starts increasing already 
at 2: ~ 20 — 25, because of the strong increase in the IGM 
temperature due to the BH emission. The predicted peak 
in STb is of the order of only a few mK, a factor of ~ 10 
smaller than in our models. The low ( 0.2) neutral fraction 
in their models is the likely cause of this discrepancy, as it 
implies a low T21 in eq. (25). 

5.3 Other X-ray feedback mechanisms 

Observations of local star-burst galaxies (Grimm, Gilfanov 
& Sunyaev 2003; Ranalli, Comastri & Setti 2003; Gilfanov, 
Grimm & Sunyaev 2004) find a correlation between star 
formation and X-ray luminosity. As was noted in Glover 
& Brand (2003), Furlanetto (2006), Pritchard & Furlan- 
etto (2007) and Santos et al. (2007), it is reasonable to ex- 
pect that also high redshift star formation is associated with 
X-ray emission, although an unknown (and possibly impor- 
tant) correction factor fx should be introduced to quantify 
the differences between the local and the primordial envi- 
ronment. 

The effects of such emission (and whether they can be 



distinguished from the ones presented in this paper) will be 
thoroughly investigated in a companion paper (Ripamonti 
et al., 2008 - in preparation). Here, we just compare the en- 
ergy input due to BHs with the one due to X-ray emission 
associated with star formation. This was done by assuming 
that the SED of such emission is a power-law with photon 
index in the 1.5 < P < 2.5 range, and that the star formation 
rate is similar to the one shown in fig. lb of Choudhury & 
Ferrara (2006). In Fig. 14 we show that in the lower redshift 
range we consider {z ^ 10) the expected energy input is only 
a fraction of the BH contribution of most of our models, and 
is comparable only to the BH model with the weakest feed- 
back; at higher redshifts the SF-associated X-ray emission 
might be more important or even dominant, but the overall 
energy input is small. 

In short, if the unknown factor fx is not much larger 
than 1, the effects of the X-ray emission associated with star 
formation should be at most comparable to the ones of the 
weakest of our BH models (such as the BVR06-f SOSl case). 



6 SUMMARY AND CONCLUSIONS 

We have examined how a population of accreting BHs might 
affect the pre-reionization Universe, looking in particular at 
the effects upon the neutral regions outside the first ionized 
'bubbles', where stellar feedback is likely small. We explored 
a number of scenarios for the growth of the cosmological 
BH density, and considered several possible SEDs. Both of 
these components are in broad agreement with observational 
constraints (e.g. the results about the X-ray background by 
Dijkstra et al. 2004). 

Our analysis started from how the energy input from 
the diffuse radiation due to the BH population might af- 
fect the temperature and ionization level of the IGM far- 
away from ionized regions where local effects are important. 
Given the Dijkstra et al. (2004) constraints, it is not sur- 
prising that BH emission in our models leads only to partial 
ionization: the main effect of BH emission is then the in- 
crease in the temperature of the IGM, which easily reaches 
levels ^ 10^ K in all the cases we have considered. 

Then, we explored a number of possible indirect conse- 
quences of the energy input: 

• CMB measurements appear unable to constrain any of 
our models, since all of them comfortably fit observational 
constraints from WMAP; 

• 21-cm observations appear extremely promising, since 
in most of the BH models the predicted STi, should be easily 
detectable with the next generation of 21-cm experiments 
(e.g. LOFAR), especially if stellar Lyman a coupling is really 
present; 

• the critical mass for halos to be able to cool, collapse 
and form stars is significantly enhanced ai z ^ 10, and in 
some of our models it becomes ~ 100 times larger than in 
the unperturbed case. This allows star formation only in ha- 
los with virial temperatures ^ 10* K, i.e. prevents (or, in 
models with weak feedback, significantly reduces) the for- 
mation of PopIII objects for z ^ 9; 

• gas depletion might occur in the models with 
intermediate-to-strong BH feedback, and for relatively low 
virialization redshifts: halos with masses between Merit and 
3—10 Merit appear to be able to form stars at their centre, 
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but their baryonic fraction is considerably lower than the 
cosmological average. 

The most relevant of our results appears to be the one 
about 21-cni observations, since it might be falsified (or 
confirmed) by forthcoming observations. To our knowledge, 
the only mechanism which should be able to heat the IGM 
outside ionized regions in a comparable way is the X-ray 
emission associated with star formation, as was proposed 
by Glover & Brand (2003), Furlanetto (2006) and Pritchard 
& FYirlanetto (2007). We leave the detailed comparison be- 
tween the two models to a future paper (Ripamonti ct al. 
2008, in preparation), where we will also investigate whether 
it is possible to distinguish between the two scenarios, e.g. 
by using the spatial power-spectrum. 

We stress that most of our conclusions strongly depend 
on the details of the reionization process, and in particular 
on the survival of neutral regions down to redshift ~ 7 — 8. 
Recent simulations (e.g. Santos et al. 2007) and observations 
(e.g. Fan et al. 2002, Kashikawa et al. 2006) hint that the 
overlapping phase lasted for a long time and suggest the ex- 
istence of patches with a significant neutral fraction even at 
z ^ 7. However, the scenario of an earlier reionization can- 
not be rejected at present: in such a case, BH signatures 
(such as the effects on the properties of 21-cm radiation) 
become difficult or impossible to detect. On the other hand, 
if the predictions of simulations are correct, the effects of BH 
emission might enhance the 21-cm contrast between nctitral 
and ionized patches, improving our capability of studying 
the 2: --^ 7 — 12 Universe, and providing important informa- 
tion on the duration and the end of the reionization phase. 



ACKNOWLEDGMENTS 

We thank R. Salvatcrra for clarifications about the CS07 re- 
sults, F. Haardt, G. Mellema, R. Thomas and M. Volonteri 
for useful discussions, and the anonymous referee for several 
suggestions about improving the manuscript. We also thank 
K. Visser for assistance in working with the HPC cluster 
of the Centre for High Performance Computing and Visu- 
alization of the University of Groningen, where most of our 
computations were carried out. ER and MM thank the In- 
stitute for Theoretical Physics of the University of Zurich, 
and the Kapteyn Astronomical Institute of the University 
of Groningen for the hospitality during the preparation of 
this paper. ER acknowledges support from the Netherlands 
Organization for Scientific Research (NWO) under project 
number 436016, and MM acknowledges support from the 
Swiss National Science Foundation, project number 200020- 
109581/1 (Computational Cosmology & Astrophysics). SZ is 
a member of the LOFAR project which is partially funded 
by the European Union, European Regional Development 
Fund, and by "Samenwerkingsverband Noord-Nederland" , 
EZ/KOMPAS. 

REFERENCES 

Abel T., Anninos P., Zhang Y., Norman M. L., 1997, NewA, 2, 
181 

Allison A. C, Dalgarno A., 1969, ApJ, 158, 423 



Anninos P., Zhang Y., Abel T., Norman M. L., 1997, NcwA, 2, 
209 

Bauer F. E., Alexander D. M., Brandt W. N., Schneider D. P., 
Treister E., Hornschemeier A. E., Garmire G. P., 2004, AJ, 
128, 2048 (B04) 

Becker R. H. et al., 2001, A J, 122, 2850 

Begelman M. C., Volonteri M., Rees M. J., 2006, MNRAS, 370, 

289 (BVR06) 
Bromm V., Locb A., 2003, ApJ, 596, 34 

dc Bruyn A. G., Zaroubi S., Koopmans L. V. E., 2007, The 
LOFAR Epoch of Reionization Project Plan, available from 
http; / /www.astro.rug.nl/LofarEoR 

de Bruyn A. G., et. al.,2008, in preparation. 

Chen X., Kamionkowski M., 2004, PhRvD, 70d, 3502 

Choudhury T. R., Ferrara A., 2006, MNRAS, 371, L55 

Ciardi B., Madau P., 2003, ApJ, 596, 1 

Ciardi B., Salvaterra R., 2007, MNRAS, in press 

(arXiv:0707.3520) (CS07) 
Grain R. A., Eke V. R., Frenk C. S., Jenkins A., McCarthy I. G., 

Navarro J. F., Pcarce F. R., 2007, MNRAS, 377, 41 
De Luca A., Molendi S., 2004, A&A, 419, 837 
Dijkstra M., Haiman Z., Loeb A., 2004, ApJ, 613, 646 
Dijkstra M., Wyithe J. S. B., Haiman Z., 2007, MNRAS, 379, 253 
Djorgovski S. G., Gastro S., Stern D., Mahabal A. A., 2001, ApJ, 

560, L5 

Eisenstein D. J., Loeb A., 1995, ApJ, 443, 11 
Fan X. ct al. 2001, AJ, 122, 2833 

Fan X., Narayanan V. K., Strauss M. A., White R. L., Becker R. 

H., Pentericci L., Rix H.-W., 2002, AJ, 123, 1247 
Fan X. et al. 2003, AJ, 125, 1649 

Fan X., Carilli C. L., Keating B., 2006, ARAA, 44, 415 
Fan X. et al. 2006, AJ, 132, 117 
Field G. B., 1958, Proc. IRE, 46, 240 

Field G. B., 1959, ApJ, 129, 525 

Furlanetto S. R., 2006, MNRAS, 371, 867 

Pritchard J. R., Furlanetto S. R., 2007, MNRAS, 376, 1680 

GaUi D., Palla F., 1998, A&A, 335, 403 

Gilfanov M., Grimm H.-J., Sunyaev R., 2004, MNRAS, 347, L57 
Glover S. G. O., Brand P. W. J. L., 2003, MNRAS, 340, 210 
Grimm H.-J., Gilfanov M., Sunyaev R., 2003, MNRAS, 339, 793 
Haardt F., Madau P., 1996, ApJ, 461, 20 
Haehnelt M. G., Rees M. J., 1993, MNRAS, 263, 168 
Hcger A., Woosley S. E., 2002, ApJ, 567, 532 
Hickox R. C, Markevitch M., 2007, ApJ, 661, L117 (HM07) 
Uicv I. T., Mellema G., Shapiro P. R., Pen U.-L., 2007, MNRAS, 
376, 534 

Islam R. R., Taylor J. E., Silk J., 2003, MNRAS, 340, 647 
Islam R. R., Taylor J. E., Silk J., 2004, MNRAS, 354, 427 
lye M. et al. 2006, Nature, 443, 186 
Kashikawa N. et al., 2006, ApJ, 648, 7 
Kogut A. et al., 2003, ApJS, 148, 161 

Koushiappas S. M., Bullock J. S., Dekel A., 2004, MNRAS, 354, 

292 

Kuhlen M., Madau P., Montgomery R., 2006, ApJ, 637, LI 
Lodato G., Natarajan P., 2006, MNRAS, 371, 1813 
Lodato G., Natarajan P., 2007, MNRAS, 377L, 64 
Loeb A., Rasio F. A., 1994, ApJ, 432, 52 

Machacek M. E., Bryan G. L., Abel T., 2003, MNRAS, 228, 273 
Madau P., Rees M. J., Volonteri M., Haardt F., Oh S. P., 2004, 

ApJ, 604, 484 
Malhotra S., Rhoads J. E., 2004, ApJ, 617, L5 
Mateo M. L., 1998, ARAA, 36, 435 
Mesinger A., Furlanetto S. R., 2007, ApJ, 669, 663 
Mesingcr A., Haiman Z., 2007, ApJ, 660, 923 
Mo H. J., Mao S., White S. D. M., 1998, MNRAS, 295, 319 
Moretti A., Gampana S., Lazzati D., Tagliaferri G., 2003, ApJ, 

588, 696 

Nusser A., 2005, MNRAS, 359, 183 



Early BH radiation - I: effects on the IGM 17 



Ostcrbrock D. E., 1989, Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei, University Science Books, Sausalito, 
CA (USA) 

Ota K. et al. 2007, arXiv0707.1561 

Pelupessy F. I., Di Matteo T., Ciardi B., 2007, ApJ, 665, 107 
Pengelly R. M., 1964, MNRAS, 127, 145 
Pritchard J. R., Purlanetto S. R., 2007, MNRAS, 376, 1680 
Ranalli P., Comastri A., Setti G., 2003, A&A, 399, 39 
Ricotti M., Ostrikcr J. P., 2004, MNRAS, 352, 547 (RO04) 
Ricotti M., Ostriker J. P., Gnedin N. Y., 2005, MNRAS, 357, 207 
Ripamonti E., Haardt F., Ferrara A., Colpi M., 2002, MNRAS, 
334, 401 

Ripamonti E., 2007, MNRAS, 376, 709 (R07) 
Ripamonti E., Mapelli M., Ferrara A., 2007, MNRAS, 375, 1399 
(RMF07) 

Salvaterra R., Haardt F., Ferrara A., 2005, MNRAS, 362L, 50 
Salvaterra R., Haardt F., Volonteri M., 2007, MNRAS, 374, 761 
Santos M. G., Amblard A., Pritchard J., Trac H., Cen R., Cooray 

A., 2007, arXiv:0708.2424 
Sazonov S. Y., Ostriker J. P., Sunyaev R. A., 2004, MNRAS, 347, 

144 

Scljak U., Zaldarriaga M., 1996, ApJ, 469, 437 

Seljak U., Sugiyama N., White M., Zaldarriaga M., 2003, Phys. 
Rev. D, 68h, 3507 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Shlosman I., Frank J., Begelman M. C., 1989, proceedings of the 
134th Symposium of the International Astronomical Union, 
held in Santa Cruz, California, August 15-19, 1988. Edited 
by Donald E. Ostcrbrock and Joseph S. Miller. International 
Astronomical Union. Symposium no. 134, Kluwcr Academic 
Publishers, Dordrecht, 462 

Shlosman I., Begelman M. C, Prank J., 1990, Nature, 345, 679 

ShuU J. M., van Steenberg M. E., 1985, ApJ, 298, 268 

Smith F.J., 1966, P&SS, 14, 929 

Spergel D. N. et al., 2003, ApJS, 148, 175 

Spergel, D. N. et al., 2007, ApJS, 170, 377 

Steidel C. C., Hunt M. P., Shapley A. E., Adelberger K. L., Pettini 
M., Dickinson M., GiavaUsco M., 2002, ApJ, 576, 653 

Tegmark M., Silk J., Rees M. J., Blancliard A., Abel T., Palla F., 
1997, ApJ, 474, 1 

Thomas R. M., Zaroubi S., 2007, MNRAS, submitted, 
(arXiv:0709.1657) 

Toomre A., 1964, ApJ, 139, 1217 

Umemura M., Loeb A., Turner E. L., 1993, ApJ, 419, 459 
Valdcs M., Ferrara A., MapeUi M., Ripamonti E., 2007, MNRAS, 

377, 245 

Volonteri M., Salvaterra R., Haardt F., 2006, MNRAS, 373, 121 
Volonteri M., Haardt F., Madau P., 2002, Ap&SS, 281, 501 
Volonteri M., Haardt F., Madau P., 2003, ApJ, 582, 559 
Volonteri M., Perna R., 2005, MNRAS, 358, 913 
Volonteri M., Rees M. J., 2005, ApJ, 633, 624 
Volonteri M., Rees M. J., 2006, ApJ, 650, 669 
White R. L., Becker R. H., Fan X., Strauss M. A., 2003, AJ, 126, 
1 

Wild J. P., 1952, ApJ, 115, 206 

Zahn O., Lidz A., McQuinn M., Dutta S., Hernquist L., Zaldar- 
riaga M., Furlanetto S. R., 2007, ApJ, 654, 12 

Zaroubi S., Thomas R. M., Sugiyama N., Silk J., 2007, MNRAS, 
375, 1269 (Z07-I-) 

Zdziarski A. A., Svensson R., 1989, ApJ, 344, 551 

Zygelman B., 2005, ApJ, 622, 1356 



